Why do bubbles in Guinness sink? 
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Stout beers show the counter-intuitive phenomena of sinking bubbles while the beer is settling. 
Previous research suggests that this phenomena is due the small size of the bubbles in these beers 
and the presence of a circulatory current, directed downwards near the side of the wall and upwards 
in the interior of the glass. The mechanism by which such a circulation is established and the 
conditions under which it will occur has not been clarified. In this paper, we demonstrate using 
simulations and experiment that the flow in a glass of stout depends on the shape of the glass. If it 
narrows downwards (as the traditional stout glass, the pint, does), the flow is directed downwards 
near the wall and upwards in the interior and sinking bubbles will be observed. If the container 
widens downwards, the flow is opposite to that described above and only rising bubbles will be seen. 
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Stout beers, such as Guinness, foam due to a com- 
bination of dissolved nitrogen and carbon dioxide^ as 
opposed to other beers which foam due to dissolved car- 
bon dioxide alone. The use of nitrogen results in a range 
of desirable characteristics of the beer, including a less 
bitter taste and a creamy long-lasting head which can be 
attributed to the low solubility of nitrogen and small size 
of the bubbles^ This small bubble size is also responsi- 
ble, at least in part, for another intriguing characteristic 
of stout beers: the phenomenon of sinking bubbles, ob- 
served while the beer is settling, i.e. between the pouring 
of the beer and the formation of the head4 

Experimental studies 5 have demonstrated that the 
phenomenon of sinking bubbles is real and not an op- 
tical illusion, while simulations^ show that the bubbles 
are driven by a downward flow, the velocity of which 
exceeds the upward velocity of the bubble due to the 
Archimedean force. The existence of such a flow near 
the wall of the glass implies that there must be an up- 
ward flow somewhere in the interior. The mechanism of 
this circulation is, however, unclear, as is the role of the 
shape of the glass. 

Understanding these types of bubbly flows is impor- 
tant for a number of applications, such as manufacturing 
champagne glasses engraved with nucleation sites £ wid- 
get and similar technologies for promoting foaming in 
stouts^, designing glasses which minimize the settling 
time of stouts and, generally, for industrial processes in- 
volving bubbly flows (e.g. bubble columns^). 

In this paper, we put forward an explanation of the 
effect of sinking bubbles in Guinness, which takes into 
account the role of the shape of the glass. First, in Sec- 
tion [n] we describe the properties of Guinness as a two- 
phase medium. In Section IIIII we present the results of 
numerical simulations for several shapes of the glass. In 
Section IIVI we explain the basic mechanism which drives 
bubbles downwards and describe a simple experiment 
which can be used to confirm our hypothesis. Finally 
we give conclusions in Section fVl 



We shall model Guinness by a liquid of density pi and 
viscosity fj>i, with randomly distributed bubbles of gas 
of density p g and viscosity p g . For a temperature of 
6°C (recommended for consumption of Guinness by its 
producer "Diageo" 12 ) and normal atmospheric pressure, 
we have 



pi = 1007 kg m" 3 
p g = 1.223 kguT 3 



LH = 2.06 x 10 _3 Pas 



Pg 



0.017 x lO^Pas 



where the former values have been measured by ourselves 
and verified against the extrapolation formula of Ref. 0. 

To check whether the bubble shapes differ from 
spheres, we introduce the Bond number 



Bo 



Pigdl 



where db is the bubbles' characteristic diameter, a is the 
surface tension of the liquid/gas interface, and g is the 
acceleration due to gravity. Assuming db = 122 |xm (as 
reported in Ref. 111! ) and a = 0.745 Nm -1 (which cor- 
responds to water/air interface), we obtain Bo ps 0.002 
- which is sufficiently small to assume that bubbles in 
Guinness are spherical. 

Note also that Guinness (as well as vast majority of 
'real' liquids) contains a lot of surfactants, which make 
the bubbles behave as rigid spheres J£ This allows one to 
estimate the characteristic bubble velocity Ub using the 
Stokes formula for a rigid sphere, 



Ub 



(Pi - Pb) gd 2 b 

ism 



3.96 mms 



Estimating the corresponding Reynolds number 
PiUbdb 



Re = 



0.24, 



confirms that the Stokes formula yields a qualitatively 
correct value for Ub- Furthermore, given that Ub is much 
smaller than the speed of sound, the gas can be treated 
as incompressible. 



2 



Finally, we introduce the void fraction /, i.e. the gas's 
share of the volume of the liquid/gas mixture. For canned 
Guinness, / « 0.05 (see Ref. Qlh . whereas for draught 
Guinness served in pubs, / rs 0.1 (according to our own 
measurements). Note, however, that, traditionally, bar- 
tenders first fill, say, 80% of the glass and wait until it 
has fully settled (i.e. all the bubbles have gone out of the 
liquid into the foamy head), after which they would fill 
the glass full. Thus, when Guinness is served to the cus- 
tomer, the void fraction can be estimated as / « 0.02, 
which is the value used in this work. 



III. NUMERICAL MODELING OF THE 
LIQUID/BUBBLE CIRCULATION 

To simulate flows in Guinness, we use the finite ele- 
ment model for bubbly flows included in the COMSOL 
Multiphysics package. The model's physical foundations 
are described in detail in Ref. In this model the bub- 
bles are assumed to be all of the same size. In view of 
the problem's axial symmetry, the axi-symmetric version 
of the model is used. 

Two geometries of the holding container were exam- 
ined (see Fig. [I]): a pint and an 'anti-pint', i.e. the pint 
turned upside-down. In both cases the initial distribution 
of bubbles was uniform, and the physical parameters of 
Guinness were as described above. 

The results of typical simulations are shown in Fig. [T] 
One can see that an elongated vortex arises near the slop- 
ing part of the pint container, resulting in a downflow of 
bubbles along the wall (see the top-left panel of Fig. A 
similar vortex also exists in the anti-pint, but it rotates in 
the opposite direction and, thus, causes an upward flow 
near the wall (see the top-right panel of Fig. [5]). 

Another important feature to be observed is the narrow 
region of low density of bubbles along the wall of the 
pint container (see the left panel of Fig. [T] and the lower- 
left panel of Fig. [5]). In the anti-pint container, in turn, 
the bubble density increases near the wall (which is not 
visible in the right panel of Fig. [TJ but can be clearly seen 
in the lower-right panel of Fig. [2]). 

We have also examined the evolution of the global void 
fraction for the pint and anti-pint, as well as the cylindri- 
cal container of the same volume. Fig. [3] shows that all 
three geometries provide more or less the same settling 
time T s (for a glass of stout, a smaller T s is generally 
regarded as an advantage). 

T s has also been used to examine the extent to which 
our results depend on the void fraction / and the bubble 
size db (these are the only parameters with 'uncertain' 
values). It has turned out that, surprisingly, the depen- 
dence of T s on / is very weak: e.g., an increase in / 
from 0.02 to 0.05 results in a decrease in T s from 43 to 
46 seconds (the settling time here was defined such that 
f(T s ) = 10 -6 ). The dependence of T s on db is much 
stronger: a decrease in db from 122 to 90 micron results 
in an increase in T s from 43 to 83 seconds. 
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FIG. 1. Numerical simulations of bubbly flows for the pint and 
anti-pint. The curves show the streamlines for the bubbles, 
the color shows the void fraction /. The snapshots displayed 
correspond to t = 4 s. Observe the region of reduced / near 
the wall of the pint (the near-wall region of increased / in the 
anti-pint is not visible in this figure, but can be observed in 
Fig. EH). 



Note also that the settling time of, approximately, 43 
seconds computed for db — 122 pm and the pint con- 
tainer does not agree with our experimental estimate of 
approximately 120 seconds. The difference between the 
two results is probably caused by the fact that all bubbles 
in our simulations were of the same size, whereas in re- 
ality they are distributed with a certain dispersion. One 
can then conjecture that smaller bubbles leave the liquid 
later than those of the median size, which would account 
for the difference between the computed and measured 
values of T s . 

In what follows, we shall argue that the circulation de- 
veloping in the flow is determined by the near- wall vari- 
ation of the bubble density (as suggested previously by 
Ref. [3 for bubble columns), and that the bubble density, 
in turn, is determined by the shape of the container. 



IV. THE MECHANISM OF THE EFFECT 

First of all, observe that, whichever way the bubbles 
move, they exert a drag force on the surrounding liquid. 
This does not mean, however, that the liquid is neces- 
sarily entrained by the motion of the bubbles. Indeed, 
if, for example, the bubbles (and, hence, the drag force) 
are distributed uniformly, all liquid particles must move 



3 



Pint Anti-pint 




q L_i i i i i i i i i I i i i i i i i i i I i i i i i i i i i I i i i i i J L_i i i i i i i i i I i i i i i i i i i I i i i i i i i i i I i i i i i l 

1 2 3 1 2 3 

r (cm) r (cm) 

FIG. 2. The half-height cross-sections of the vertical velocity u and the void fraction / for the pint and anti-pint geometries 
(these graphs correspond to the (r, z) diagrams shown in Fig. 2). The dotted lines in the upper panels separate the regions of 
upward/downward flow. 




t (s) 

FIG. 3. The global void fraction / (i.e. the proportion of 
the gas in the container) versus time t, for the three cases: 
the pint (p), anti-pint (a) — both illustrated in Fig. [2] — and a 
cylinder of the same volume (c). 



the same way - which effectively means that they cannot 
move at all due to the liquid's incompressibility and the 
fact that the container has a bottom. In this case, the 
drag force is compensated by a pressure gradient exerted 
in the fluid. 

Let us now assume that there is a region of low bubble 
density near the container's wall (as there indeed is in 
the pint container). In this case, the density of the drag 
force near the container's axis is larger than that near the 
wall - which creates an imbalance and, thus, gives rise 
to a circulation: near the axis, the liquid flows upwards 
and, near the wall, downwards. Then, if the velocity of 



the downward flow is larger than the relative velocity 
of the bubbles, the bubbles will be observed to sink. A 
similar argument indicates that a near-wall region with 
higher bubble density gives rise to an upward flow (i.e. 
exactly what our simulations show for the anti-pint con- 
tainer) . 

It still remains to identify the mechanism reducing the 
bubble density near the wall for the pint geometry and 
increasing it for the anti-pint one. Potentially, there are 
two such mechanisms. The most obvious one is based 
on the 'lift force', generated by the flow around a sphere 
moving along a rigid boundary: this force pushes the 
sphere away from the boundary^. In the limit of small 
Reynolds number, however, the lift force is weak, and 
estimates show that the resulting reduction of the bubble 
density near the wall is negligible. Another shortcoming 
of this mechanism is that it docs not seem to distinguish 
between the pint and anti-pint geometries. 

To explore the effect of the geometry, assume that the 
container is not cylindrical, but narrows slightly towards 
its bottom (as the pint does). Then, even if the bubbles 
were initially distributed uniformly, their upward motion 
immediately creates a bubble-free zone along the wall 
(see Fig. 2]). On the other hand, in a container that 
widens towards its bottom (as the anti-pint does), the 
initially upward motion of bubbles increases the near- 
wall bubble density. We believe that this simple kine- 
matic effect is responsible for the circulation observed in 
Guinness. 

This effect, although not previously discussed in the 
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FIG. 4. The evolution of bubbles near the wall (a schematic), 
for (a) container narrowing downwards (the bubbles move 
away from the wall); (b) container widening downwards (the 
bubbles move towards the wall and, eventually, accumulate 
there) . 



context of Guinness, is well known in sedimentation the- 
ory as the Boycott effect liLiS It was first observed in test 
tubes containing red blood cells when it was discovered 
that sedimentation times could be significantly reduced 
by inclining the test tubes. 



Finally, our conclusions can be readily verified exper- 
imentally. If Guinness is poured into a tall cylindrical 
container (such as, for example, a laboratory measuring 
cylinder) and if the container is tilted, bubbles will be 
observed to move upwards near its upper surface and 
downwards near its lower surface - in precise agreement 
with the mechanism proposed (see supplementary mate- 
rial for a video of this experiment). 



V. CONCLUSIONS 

The sinking bubbles of Guinness and other stout beers 
have intrigued beer drinking physicists and their students 
for some time. In this paper we describe the role the 
shape of the Guinness pint glass plays in promoting the 
circulatory flow responsible for the sinking bubbles. In 
doing so we complete the explanation of this phenom- 
ena, building on previous experimental and simulation 
work. This allows us to understand the physics underly- 
ing the shape of the Guinness pint glass. It also raises the 
intriguing question: is this shape the most efficient pos- 
sible or could the settling time be significantly reduced 
by some other, possibly non-axisymmetric, shape of pint 
glass? 
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